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Abstract. The hydrodynamics of a variant of classical 
Bondi-Hoyle-Lyttleton accretion is investigated: a totally 
absorbing sphere moves at various Mach numbers (3 and 
10) relative to a medium, which is taken to be an ideal gas 
having a density gradient (of 3%, 20% or 100% over one 
accretion radius) perpendicular to the relative motion. I 
examine the influence of the Mach number, the adiabatic 
index, and the strength of the gradient upon the physical 
behaviour of the flow and the accretion rates of the angular 
momentum in particular. The hydrodynamics is modeled 
by the "Piecewise Parabolic Method" (PPM). The resolu- 
tion in the vicinity of the accretor is increased by multiply 
nesting several grids around the sphere. 

Similarly to the 3D models published previously, both 
with velocity gradients and without, the models with a 
density gradient presented here exhibit non-stationary 
flow patterns, although the Mach cone remains fairly sta- 
ble. The accretion rates of mass, linear and angular mo- 
menta do not fluctuate as strongly as published previ- 
ously for 2D models. No obvious trend of the dependency 
of mass accretion rate fluctuations on the density gradient 
can be discerned. The average specific angular momentum 
accreted is roughly between zero and 70% of the total an- 
gular momentum available in the accretion cylinder in the 
cases where the average is prograde. Due to the large fluc- 
tuations during accretion, the average angular momentum 
of some models is retrograde by up to 25%. The magnitude 
is always smaller than the value of a vortex with Kepler 
velocity around the surface of the accretor. 

The models with small density gradients initially dis- 
play a transient quasi-stable accretion phase in which the 
specific angular momentum accreted is within 10% of the 
total angular momentum available in the accretion cylin- 
der. Later, when the flow becomes unstable, the average 
decreases. I conclude that for accretion from a medium 
with both density and/or velocity gradients, most of the 
angular momentum that is available in the accretion cylin- 
der is accreted together with mass. Small gradients hardly 
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influence the average accretion rates as compared to accre- 
tion from a homogeneous medium, while very large ones 
succeed to dominate and form an accretion flow in which 
the sense of rotation is not inverted. 

Key words: Accretion, accretion disks - Hydrodynam- 
ics - Instabilities - Shock waves - Methods: numerical - 
Binaries: close 



1. Introduction 

The simplicity of the classic Bondi-Hoyle-Lyttleton (BHL) 
accretion model makes its use attractive in order to esti- 
mate roughly accretion rates and drag forces applicable 
in many different astrophysical contexts. Various aspects 
of the BHL flow have repeatedly been investigated in the 
past by many authors. In the BHL scenario a totally ab- 
sorbing sphere of mass M moves with velocity Voa rel- 
ative to a surrounding homogeneous medium of density 
Poo and sound speed Coq. In this second instalment, I ex- 
tend the investigations started in Ruffcrt (1997, hence- 
forth Rl) to include density gradients of the surrounding 
medium. Usually, the accretion rates of various quanti- 
ties, like mass, angular momentum, etc., including drag 
forces, are of interest as well as the bulk properties of 
the flow, (e.g. distribution of matter and velocity, stabil- 
ity, etc.). The results pertaining to total accretion rates 
tend to agree well qualitatively (to within factors of two, 
and ignoring the instabilities of the flow) with the original 
calculations of Bondi, Hoyle and Lyttleton (Hoyle & Lyt- 
tleton 1939, 1940a, 1940b, 1940c; Bondi & Hoyle 1944). 
However, the question of whether and how much angular 
momentum is accreted together with mass from an in- 
homogeneous medium has remained largely unanswered, 
although Rl has attempted a first answer. It has already 
been summarised in the introduction of Rl, that the strict 
application of the BHL recipe to inhomogeneous media, 
including some small constant gradient in the density or 
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Table 1. Parameters and some computed quantities for all models. Moo = Moo(y = 0) is the Mach number of the unperturbed 
flow, e p the parameter specifying the magnitude of the density gradient, 7 the ratio of specific heats, tf the total time of the run 

(units: Ra/coo), M the integral average of the mass accretion rate, S one standard deviation around the mean M of the mass 

accretion rate fluctuations, M the maximum mass accretion rate, Mbh is defined in Eq. (3) of Ruffert & Arnett (1994), Z x , ly, 
l z , are the averages of specific angular momentum components together with their respective standard deviations u x , cr y , a z , s 
is the entropy (Eq. (4) in Ruffert & Arnett 1994), the radius of the accretor always is R± — 0.02J?a, the number of grid nesting 
depth levels g = 9, the size of one zone on the finest grid S — l/256i?A, the softening parameter for the potential of the accretor 
(sec Ruffert 1994) e = 3 zones, the number N of zones per grid dimension is 32, and the size of the largest grid is L — 32Ra- 
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the velocity distribution, yields that the accreted matter 
has zero angular momentum by construction (Davies & 
Pringle, 1980). 

We recall that the largest radius i? a from which matter 
is still accreted by the BHL-procedure turns out to be the 
so-called Hoyle-Lyttleton accretion radius 



where G is the gravitational constant. The mass accretion 
rate follows to be 

M HL = nRlpooVoo ■ (2) 

I will refer to the volume upstream of the accretor from 
which matter is accreted as accretion cylinder. One can 
additionally calculate (Dodd & McCrea 1952; Illarionov 
& Sunyaev 1975; Shapiro & Lightman 1976; Wang 1981) 
how much angular momentum is present in the accretion 
cylinder for a non-axisymmetric flow which has a gradient 
in its density or velocity distribution perpendicular to the 
mean velocity direction. Then, assuming no redistribution 
of angular momentum, the amount accreted is equal to 
(or at least is a large fraction of) the angular momentum 
present in the accretion cylinder. 

The uncertainty about how much angular momentum 
can actually be accreted in a BHL flow stems from these 
two opposing views involving either a large or a very small 



fraction of what is present in the accretion cylinder. In 
the first paper Rl, I showed that the answer is not clear 
cut, but depends on the initial and boundary conditions. 
Roughly 7% to 70% of the total angular momentum avail- 
able in the accretion cylinder is accreted. 

In this second paper I would like to compare the 
accretion rates of the angular momentum of numeri- 
cally modeled accretion flows with density gradients to 
the previous results of accretion with velocity gradients 
(Rl). Although several investigations of two-dimensional 
flows with gradients exist (Anzer et al. 1987; Fryxell & 
Taam 1988; Taam & Fryxell 1989; Ho et al. 1989; Benen- 
sohn et al 1997, Shima et al, 1998), i/iree-dimensional sim- 
ulations are scarce due to their inherently high computa- 
tional load. Livio et al. (1986) first attempted a three- 
dimensional model including gradients, but due to their 
low numerical resolution the results were only tentative. 
In the models of Ishii et al. (1993) the accretor was only 
coarsely resolved, while the results of Boffin (1991) and 
Sawada et al. (1989) are only indicative, because due to 
the numerical procedure the flows remained stable (too 
few SPH particles in Boffin 1991 and local time stepping 
in Sawada et al. 1989 which was described to be appropri- 
ate only for stationary flows). Also Sawada et al. (1989) 
only investigated velocity gradients. 

In order to be able to compare the new results to pre- 
vious models of Rl , I will use the same values for the gra- 



M. Ruffert: Non-axisymmetric wind-accretion simulations 



3 



2.1. Numerical Procedure 
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Fig. 1. The coefficient for the mass accretion rate / m (dotted), 
defined by Eq. (^), and for the specific angular momentum /j 
(solid), defined by Eq. as a function of e p . The thin curves 
show the values of / for a simple linear relation between v xoo 
and E p , while the bold curves apply to a relation including the 
"tanh" -term as given in Eq. (Q). The two vertical straight lines 
indicate the gradients that were used in the numerical models 
(cf. Table 0), e p = 0.03 and e p = 0.2. Compare this figure with 
Fig. 1 in paper Rf and note the differences. 



dients as in Rl. However, as will be mentioned in connec- 
tion with Eq. ^| (which had already been presented in Rl), 
the expected angular momentum accretion is reduced by 
a factor of six, when changing from velocity to density 
gradients. Thus the clear separation between angular mo- 
mentum from the bulk flow one from the unstable, fluc- 
tuating flow blurs. On the other hand arbitrarily large 
specific angular momenta cannot be accreted because of 
the ang. mom. barrier felt by matter spiralling into the 
accretor. A simulation with a very large density gradient 
was done, too, mainly to be able to compare to a previous 
2D model (Fryxell & Taam, 1988). 

In Sect. |^ I give an only short summary of the numer- 
ical procedure used. Sections [| and ^ present the results, 
which I summarise in Sect. [| 



2. Numerical Procedure and Initial Conditions 

Since the numerical procedures and initial conditions are 
mostly identical to what has already been described and 
used in previous papers (cf. Ruffert 1997, Rl, and refer- 
ences therein) I will refrain from repeating every detail, 
but only give a brief summary. 



The distribution of matter is discretised on multiply 
nested equidistant Cartesian grids (e.g. Berger & Colella, 
1989) with zone size S and is evolved using the "Piece- 
wise Parabolic Method" (PPM) of Colella & Woodward 
(1984). The equation of state is that of a perfect gas with a 
specific heat ratio 7 (see Table |l|) . The model of the maxi- 
mally accreting, vacuum sphere in a softened gravitational 
potential is summarised in Ruffert & Arnett (1994) and 
Ruffert & Anzer (1995). 

A gravitating, totally absorbing "sphere" moves rela- 
tive to a medium that has a distribution of density and 
velocity far upstream (at x — > —00) given by 
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with the redefined accretion radius (note slight difference 
to Eq. 0) 



Ra = 



2GM 



(5) 



In this paper I only investigate models with gradients 
of the density distribution, thus for all models I set e v = 0. 
The values of s p can be found in Table [l| In order to keep 
the pressure constant throughout the upstream boundary, 
I changed the internal energy or equivalently the sound 
speed c accordingly (c 2 = jp/p)- Thus the Mach number 
.Moo (?/) = v/coo(y) of the incoming flow varies too, since 
the velocity is kept constant. The Mach numbers given 
in the Table [l] refer to the values along the axis (y=0) 
upstream of the accretor. 

The function "tanh" is introduced in Eq. (g|) and (|J) 
to serve as a cutoff at large distances y for large gradients 
e. The units I use in this paper are (1) the on-axis sound 
speed Coo(y = 0) as velocity unit; (2) the accretion radius 
(Eq. (H)) as unit of length, and (3) the on-axis density po 
. Thus the unit of time is -Ra/coo- 

I'll assume both e v 7^ and e p 7^ for the following es- 
timates. Accreting all mass within the accretion cylinder, 
and taking the matter to have the density and velocity 
distributions as given by Eqs. (§) and (|) one obtains the 
mass accretion rate to lowest order in e v and e p to be 



M = ttR a pqv q 



(6) 



an equation very similar to Eq. (g). Further assuming 
that all angular momentum within the deformed accretion 
cylinder is accreted too, the specific angular momentum of 
the accreted matter follows to be (Ruffert & Anzer 1994; 
Shapiro & Lightman 1976; again to lowest order in e v ) 



1 



(6e v + e p ) v Ra 



(7) 
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Fig. 2. Contour plots showing snapshots of the density to- 
gether with the flow pattern for all models with an adiabatic 
index of 5/3. The contour lines are spaced logarithmically in 
intervals of 0.1 dex. The bold contour levels are sometimes la- 
beled with their respective values (0.0, 0.2, and 0.4). Darker 
shades of gray indicate higher densities. The dashed contour 
delimits supersonic from subsonic regions. The time of the 
snapshot together with the velocity scale is given in the leg- 
end in the upper right hand corner of each panel. 




Fig. 3. The accretion rates of several quantities are plotted as a function of time for model MV (A^oo^l.4, gradient of 3%, 
adiabatic index of 7 = 5/3). The left panel contains the mass and angular momentum accretion rates, the right panel the 
specific angular momentum of the matter that is accreted. In the left panel, the straight horizontal lines show the analytical 
mass accretion rates: the top solid line is the Bondi-Hoyle approximation formula (Eq. (3) in Ruffert 1994; Bondi 1952) and the 
lower one half that value. The upper solid bold curve represents the numerically calculated mass accretion rate. The lower three 
curves of the left panel trace the x (dotted), y (thin solid) and z (bold solid) component of the angular momentum accretion 
rate. The same components apply to the right panel. The horizontal line in the right panel shows the specific angular momentum 
value as given by Eq. (pi). 



For positive e p the density is lower on the positive side of 
the y-axis, then the vortex formed around the accretor is 
in the anticlockwise direction, i.e. the angular momentum 
component in z-direction is positive. 

The values of the specific angular momentum obtained 
from the numerical simulations will be compared to the 
values that follow from this Eq. (Q) to conclude which of 
the above mentioned views — low or high specific angular 
momentum of the accreted material — is most more ap- 
propriate. As was stated further above, Rl found that a 
sizable amount (between 7% and 70%) is accreted, when 
velocity gradients are present. I will implicitly assume the 
component j z when discussing properties like fluctuations, 
magnitudes, etc. From the symmetry of the boundary con- 
ditions the average of the x and y components of the angu- 
lar momentum should be zero, although their fluctuations 
can be quite large. Taking the numerically obtained ac- 
cretion rates of mass M(t) and angular momentum J(t), 
I calculate the instantaneous specific angular momentum 
j(t) = J(t)/M(t), the temporal mean I of which is listed 
in Table fjj too. 

Apart from serving as cut-off, the tanh-dependencies in 
Eqs. and have a gradient that is slightly less steep 
than simply linear. Thus less specific angular momentum 
is present at a given radius from the accretor and smaller 
values in the magnitude of the accreted quantity result. 

One can numerically approximate the integrals of the 
mass flux and angular momentum over the cross section 



of the accretion cylinder, to obtain the coefficients / in 
the relations Eq. (||) and Eq. (Q) : 

M = f m (ep)TTR 2 A p vo , (8) 

Jz = fi(e P )e p voRA • (9) 

Here, I will only consider the effect of a density gradient. 
The unitless functions / are a function of e p and the func- 
tional relation of /foo[e p ], i.e. whether depends purely 
linearly on e p or as in Eq. (||) via the "tanh"-term. Fig- 
ure [l] shows the values of the functions / for the mass and 
specific angular momentum and for both the linear and 
"tanh" case. The only minimal deviation of the "tanh" 
curves from the linear ones indicates that the tanh-cutoff 
hardly acts within the accretion cylinder. Since f m rj 1 
and is practically constant for s p ^ 1, Eq. (||) is a good 
approximation in this range. If the prescription is cor- 
rect that everything in the accretion cylinder is accreted, 
no difference in the accretion rate should be observable 
in the models with differing magnitude of gradient e p 
(cf. Table 0). The same constancy applies to the accre- 
tion of specific angular momentum: its coefficient /j re- 
mains relatively constant /j rj 0.25 in the range e p ^ 0.4, 
which includes both e p for which models were simulated. 
So although I used the same magnitudes for the gradients, 
the effects I expect in the accretion rates are markedly 
different from what was presented in paper Rl. Only for 
model VS with e p — 1.0, does /j deviate appreciably from 
0.25: it is /j w 0.16. 
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Fig. 4. The accretion rates of several quantities are plotted as a function of time for the moderately supersonic (.Moo=3) 
models MS and NS with an index 7 = 5/3. The top panels contain the mass and angular momentum accretion rates, the 
bottom panels the specific angular momentum of the matter that is accreted. In the top panels, the straight horizontal lines 
show the analytical mass accretion rates: half the value of the Bondi-Hoyle approximation formula (Eq. (3) in Ruffert 1994; 
Bondi 1952). The upper solid bold curve represents the numerically calculated mass accretion rate. The lower three curves of 
the top panels trace the x (dotted), y (thin solid) and z (bold solid) component of the angular momentum accretion rate. The 
same components apply to the bottom panels. The horizontal line in the bottom panels show the specific angular momentum 
value as given by Eq. 



2.2. Models 

The combination of parameters that I varied, together 
with some results are summarised in Table [IJ. The first 
letter in the model designation indicates the strength of 
the gradient: M, P, and T have e p = 0.03, while N, Q, and 
U have e p = 0.2 and V is e p = 1.0. The second letter spec- 
ifies the relative wind flow speeds, F (fast), S (slow) and 
V stand for Mach numbers of 10, 3 and 1.4, respectively. 
I basically simulated models with all possible combina- 
tions of the two higher flow speeds (Mach numbers of 3 
and 10), the two gradients (3% and 20%) and varying the 
adiabatic index between 5/3, 4/3 and 1.01. Model VS with 



the largest gradient of 100% facilitates a comparison to a 
previous two-dimensional simulation by Fryxell & Taam 
(1988). The grids are nested to a depth g = 9 such that the 
radius of the accretor i?* = 0.02i?A spans several zones on 
the finest grid. 

As far as computer resources permitted, I aimed at 
evolving the models for at least as long as it takes the flow 
to move from the boundary to the position of the accretor 
which is at the centre (crossing time scale). This time is 
given by L/2A / l o and ranges from about 1 to about 10 
time units. The actual time if that the model is run can 




Fig. 5. The accretion rates of several quantities are plotted as a function of time for the highly supersonic (.Moo = 10) models MF 
and NF with an adiabatic index of 7 = 5/3. The top panels contain the mass and angular momentum accretion rates, the bottom 
panels the specific angular momentum of the matter that is accreted. In the top panels, the straight horizontal lines show the 
analytical mass accretion rates: half the value of the Bondi-Hoyle approximation formula (Eq. (3) in Ruffert 1994; Bondi 1952). 
The upper solid bold curve represents the numerically calculated mass accretion rate. The lower three curves of the top panels 
trace the x (dotted), y (thin solid) and z (bold solid) component of the angular momentum accretion rate. The same components 
apply to the bottom panels; the horizontal line shows the specific angular momentum value as given by Eq. (pi). 



be found in Table [j], as well as the parameters of the grids 
(L, g, etc.). 

When modeling a BHL flow with a density gradient, 
one has to pay attention to the fact that matter parcels 
with possibly very different densities (which initially are 
separated upstream) will be focussed to find themselves 
close to each other along the accretion axis (U. Anzer, 
personal communication). This might render the flow ad- 
ditionally unstable. Although true in principle, this prob- 
lem does not affect, in practice, the models I will present: 
a closer inspection of Eq. (|^) reveals that the density jump 
from one end of the accretion radius to the other is only a 
factor 1.5 for the case with large gradient (e = 0.2). This 



does not seem to have an additional influence as compared 
to models with constant density presented in paper Rl. 

The calculations are performed on a Cray-YMP 4/64 
and a Cray J90 8/512. They need about 12-16 MWords 
of main memory and take approximately 160 CPU-hours 
per simulated time unit. 

3. Dynamics and accretion rates 

3.1. Results of models with "/=5/3 and 4/3 

I will describe the results of the models for which a ratio 
of specific heats of 7=5/3 was chosen together with those 
models of 7=4/3 because the evolution is very similar. The 
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Fig. 6. Contour plots showing snapshots of the density together with the flow pattern in a plane containing the centre of the 
accretor for all models with an adiabatic index of 4/3. The contour lines are spaced logarithmically in intervals of 0.1 dex. The 
bold contour levels are labeled with their respective values (0.0 or 1.0). Darker shades of gray indicate higher densities. The 
dashed contour delimits supersonic from subsonic regions. The time of the snapshot together with the velocity scale is given in 
the legend in the upper right hand corner of each panel. 



only exception is model MV: its slow relative bulk velocity 
of Mach 1.4 produces a significantly more stable flow. 

Figs. H and ^| show snapshots of the flow velocities and 
density distribution in the x-y-plane containing the accre- 
tor. The velocity pattern as well as the density contours 
within the shock cone indicate a strongly unstable flow. 
Also the shock cone itself has many bumps and kinks. 

Note how the density contours strongly bend over up- 
stream of the shock (at x ~ —0.8) for models N and Q 



indicating the large gradient as compared to models M 
and P, where only the contour of p = 1 is seen to be 
detached from the shock. The higher densities on the 
lower side of y (where it is negative) seem to produce 
an asymmetric shock cone: the lower part of the cone 
subtends a smaller angle with the y — 0-axis than the 
upper half. However, the quickly varying accretion flow 
pattern tends to impinge on the shock front and dislodge 
it. This asymetry has been observed and commented on in 
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Fig. 7. The accretion rates of several quantities are plotted as a function of time for the moderately supersonic (.Moo =3) 
models PS and QS with an adiabatic index of 7 = 4/3. The top panels contain the mass and angular momentum accretion rates, 
the bottom panels the specific angular momentum of the matter that is accreted. In the top panels, the straight horizontal lines 
show the analytical mass accretion rates: dotted is the Hoyle-Lyttleton rate (Eq. (1) in Ruffert 1994), solid is the Bondi-Hoyle 
approximation formula (Eq. (3) in Ruffert 1994; Bondi 1952) and half that value. The upper solid bold curve represents the 
numerically calculated mass accretion rate. The lower three curves of the top panels trace the x (dotted), y (thin solid) and 
z (bold solid) component of the angular momentum accretion rate. The same components apply to the bottom panels; the 
horizontal line shows the specific angular momentum value as given by Eq. (m). 



Ishii et al. (1993) and Soker & Livio (1984). In the models 
with smaller 7 (models P and Q) the shock cone tends to 
be narrower and its minimum distance from the accretor 
smaller than for the models M and N. 

In contrast to what has been said, model MV with 
a slow (but still supersonic relative velocity), shows a 
very regular flow pattern: the stagnation point is about 
0.3 Ra downstream from the accretor. Matter that comes 
within this point gets accreted while matter that stays 
outside just passes the accretor. The mass accretion rate 
of model MV (cf. Fig. ||) rises slowly and nearly monoton- 
ically to saturate close to the Bondi-Hoyle formula value. 



Only towards the end is there an indication that the flow 
might become unstable for this model too. The compo- 
nent of interest (z) for the specific angular momentum 
saturates at about 75% of the value estimated analyti- 
cally by Eq. (0). This is probably due to the fact that for 
bulk velocities close to a Mach number of unity the Hoyle- 
Lyttleton approximation brakes down because a certain 
fraction of mass is accreted practically spherically sym- 
metrically, i.e. as Bondi tried to describe it with the ap- 
proximative formula (Bondi, 1952). The other two compo- 
nents (x and y) fluctuate around zero (which is given by 
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Fig. 8. The accretion rates of several quantities are plotted as a function of time for the highly supersonic (.Moo=10) models PF 
and QF with an adiabatic index of 7 = 4/3. The top panels contain the mass and angular momentum accretion rates, the 
bottom panels the specific angular momentum of the matter that is accreted. In the top panels, the straight horizontal lines 
show the analytical mass accretion rates: dotted is the Hoyle-Lyttleton rate (Eq. (1) in Ruffert 1994), solid is the Bondi-Hoyle 
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z (bold solid) component of the angular momentum accretion rate. The same components apply to the bottom panels; the 
horizontal line shows the specific angular momentum value as given by Eq. (ph. 



the symmetry of the boundary conditions), but indicate 
at a very early stage that the flow is mildly unstable. 

The other 8 models (MS, MF, NS, NF, PS, PF, QS, 
QF) show very strong fluctuations of the accretion rates of 
mass and all angular momentum components (Figs. ||, || 
0, and ||) . A variation of factors of two is not uncommon, 
so the averages stated in Table [j] should be used with care 
and where possible the standard deviations taken into ac- 
count. I include the accretion rate plots for all models in 
order to facilitate the judgement of how representative the 
average values are for the whole temporal evolution. 



A few trends can be discerned. All four models MS, 
MF, PS, PF, i.e. the ones with small e p = 0.03 dis- 
play a fairly quiet initial transient phase. The x and y- 
components of the angular momentum fluctuate mildly 
around zero and the z-component reaches fairly precisely 
(to within 10%) the analytic estimate Eq. (Q). This shows 
the applicability of the analytic result only to the initial 
quasi-stationary state. The models with large e p = 0.2 
do not have this quiet phase but become chaotic much 
quicker. This decreased stability also produces lower ac- 
cretion rates for mass as well as lower specific angular 
momenta. The initial quiet transient phase has already 
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Fig. 9. Left panels: The accretion rates of several quantities 
are plotted as a function of time for the model with largest 
gradient (e p = 1.0). The top panel contains the mass and 
angular momentum accretion rates, the bottom panels the 
specific angular momentum of the matter that is accreted. In 
the top panel, the straight horizontal line shows the analyt- 
ical half the Bondi-Hoyle value from the approximation for- 
mula (Eq. (3) in Ruffert 1994; Bondi 1952). The upper solid 
bold curve represents the numerically calculated mass accre- 
tion rate. The lower three curves of the top panel trace the 
x (dotted), y (thin solid) and z (bold solid) component of the 
angular momentum accretion rate. The same components ap- 
ply to the bottom panel. Right panel: Comparison of the 
ratio of specific angular momenta accreted for models with 
density gradients (MS, MF, PS, NS, NF) to models with ve- 
locity gradients (IS, JS, SS, KS, LS taken from paper Rl), 
keeping all other parameters equal. 



been seen in Rl: compare the present models to the ones 
from Rl, IS, JS, KS, LS in Figs. 4, 6, 8, and 9 in Rl, 
respectively. 

Note that model PF is one for which the simulation 
was run fairly long compared to the timescale of fluctua- 
tions. This increases the confidence that the average mass 
accretion rate is a representative value and not a random 
one of a transient state. The angular momenta, however, 
still do not display the marginal positive shift that the 
z-component should have compared to the other two. 

No significant difference is observed when comparing 
two models that only differ in Mach number. However the 
mass accretion rates are significantly larger in the mod- 
els with smaller 7, again ceteris paribus. The same ap- 
plied to the velocity gradient models of Rl , cf . model IS 
and SS in Figs. 4 and 11 respectively. This dependence of 
the accretion rate on the adiabatic index 7 is well known 
for stationary flows from analytic calculations, e.g. Bondi 



(1952) for spherically symmetric flows and Sect. 5 and 6 
and Figs. 9 and 10 in Foglizzo & Ruffert (1997). 

3.2. Results of models with j—l.Ol 

The first main obvious difference between the nearly 
isothermal models and the more adiabatic ones is that 
the shock cone is attached to the accretor, as can be seen 
in Fig. [l(| The pressure around the accretor in this case 
is not sufficient to push away and support the shock cone. 
Not even the large density gradient is able to dislodge the 
shock from touching the surface of the accretor. 

Both models with high Mach number (TF and UF) 
hardly show any activity of unstable flow within the shock 
cone, contrary to the moderately supersonic cases (TS and 
US). Again (as in model MV) a clear and stable stagnation 
point is present downstream of the accretor for these quiet 
models. It is a common feature of practically all simula- 
tions (including the ones by other authors) that when the 
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Fig. 10. Contour plots showing snapshots of the density together with the flow pattern in a plane containing the centre of the 
accretor for all models with an adiabatic index of 1.01 The contour lines are spaced logarithmically in intervals of 0.1 dex. The 
bold contour levels are labeled with their respective values (0.0 or 1.0). Darker shades of gray indicate higher densities. The 
dashed contour delimits supersonic from subsonic regions. The time of the snapshot together with the velocity scale is given in 
the legend in the upper right hand corner of each panel. 



shock cone is attached to the accretor no (or hardly any) 
instability is observed. Whether or not the shock cone is 
attached depends on many physical attributes of the mod- 
els (e.g. stiffness of the equation of state, etc.) and numeri- 
cal parameters (not least resolution). However, when these 
conditions collude to produce an attached shock, invari- 
ably the flow remains stable. 

The very active flow of the slower models (TS and 
US) reflects itself in a higher variability of the mass ac- 



cretion rate as compared to models TF and UF (Figs. [Tl| 
andjjJJ). However the fluctuation of the mass accretion 
rates of all these 7 = 1.01 models is much smaller than 
the fluctuations shown by the more adiabatic models de- 
scribed further above. The average mass accretion rates 
of the slow models slightly exceed the values predicted by 
the Hoyle-Lyttleton theory, while the faster models only 
reach 80% of Mbh- The z-component of the specific an- 
gular momentum in the models with large gradient (US, 
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Fig. 11. The accretion rates of several quantities are plotted as a function of time for the moderately supersonic (.Moo =3) 
models TS and US with an adiabatic index of 7 = 1.01. The top panels contain the mass and angular momentum accretion rates, 
the bottom panels the specific angular momentum of the matter that is accreted. In the top panels, the straight horizontal lines 
show the analytical mass accretion rates: dotted is the Hoyle-Lyttleton rate (Eq. (1) in Ruffert 1994), solid is the Bondi-Hoyle 
approximation formula (Eq. (3) in Ruffert 1994; Bondi 1952) and half that value. The upper solid bold curve represents the 
numerically calculated mass accretion rate. The lower three curves of the top panels trace the x (dotted), y (thin solid) and 
z (bold solid) component of the angular momentum accretion rate. The same components apply to the bottom panels; the 
horizontal line shows the specific angular momentum value as given by Eq. (m). 



UF) practically never reach the analytical estimate Eq. (^) 
while the models with small gradient (TS, TF) fluctuate 
strongly about this value, exceeding it by a factor of three 
or even reversing sign. 

These models can be compared to models without 
gradients (Ruffert, 1996) but with the same remaining 
parameters: models TF and UF should be compared to 
model HS in Fig. 7e, while models TS and US can be com- 
pared to model GS in Fig. 5e of Ruffert (1996). Both the 
different behaviour of the mass accretion rates as well as 
the amplitudes of fluctuation of the angular momentum 
accretion rates are comparable between the models, indi- 



cating that the presence of a density gradient does not sig- 
nificantly alter the accretion properties of such a strongly 
unstable flow. Of course, the mean of the z-component is 
not zero in flows with gradients. Fig. [ll] can also directly 
be compared with Fig. 11 from Ishii et al. (1993). Both the 
mass and angular momentum accretion rates show very 
similar magnitudes and evolution. 



3.3. Results of model VS 

The density gradient in the previous set of models of 



0.03 and e c 



0.2 was chosen in order to facili- 
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Fig. 12. The accretion rates of several quantities are plotted as a function of time for the highly supersonic (A4oo = 10) models TF 
and UF with an adiabatic index of 7 = 1.01. The top panels contain the mass and angular momentum accretion rates, the 
bottom panels the specific angular momentum of the matter that is accreted. In the top panels, the straight horizontal lines 
show the analytical mass accretion rates: dotted is the Hoyle-Lyttleton rate (Eq. (1) in Ruffert 1994), solid is the Bondi-Hoyle 
approximation formula (Eq. (3) in Ruffert 1994; Bondi 1952) and half that value. The upper solid bold curve represents the 
numerically calculated mass accretion rate. The lower three curves of the top panels trace the x (dotted), y (thin solid) and 
z (bold solid) component of the angular momentum accretion rate. The same components apply to the bottom panels. The 
horizontal line in the bottom panels show the specific angular momentum value as given by Eq. (M). 



tate the direct comparison between the results presented 
in this paper and the ones with velocity gradients shown 
in Rl. However, Eq. fj] then predicts that the specific an- 
gular momentum accreted will be six times smaller for the 
density gradients as compared to accretion with velocity 
gradients. That this is true becomes clear when compar- 
ing the plots for models MS, MF, PS, NS and NF with 
equivalent models from paper Rl: IS, JS, SS, KS and LS, 
respectively: In the models presented in this paper, the 
z-component of the angular momentum, which is the one 
influenced by the gradient, hardly rises above the ran- 
dom fluctuations of the other two components. In order 



to check the correct separation of the effect of the density 
gradient from the unstable nature of the flow, a model 
with larger gradient is helpful and will be presented in 
this subsection. 

The particular value of e p = 1 is suggested because 
for this case, Taam & Fryxell (1988) have found that a 
quasi-steady disk forms which does not change its sense 
of direction. Fig. || (left panels) shows the temporal evolu- 
tion of the mass and angular momentum accretion for this 
large density model. The mass accretion rate remains very 
low over the whole simulated time on average being less 
than half the value of model QS (which is similar to VS in 
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all parameters except e p ). Note that these average values 
listed in Table [l] are normalised to Eq. and not to Eq. [|. 
Thus a big part of the reduction of specific angular mo- 
mentum between model QS (0.55) and model VS (0.23) is 
probably due to the "tanh"-term represented by the fac- 
tor fj (Fig. 0): fj w 0.25 for model QS while fj « 0.16 for 
model VS. 

Also note that model VS does not display as drastic 
a reduction in mass accretion rate as shown in Fig. 22 of 
Taam & Fryxell (1988) [note the different choice of time 
units between this paper and the one in Taam & Fryxell, 
1988]. The reason probably being that in their 2D models 
accretion gets effectively shut off as soon as a stable disk 
structure forms, while in my 3D simulations accretion can 
still proceed practically unimpeded via the poles. 

The lower left panel of Fig. || shows the angular mo- 
mentum accretion. In this large gradient model, the z- 
component clearly dominates compared to the other two 
components, i.e. the fluctuating flow cannot compete with 
the angular momentum available in the bulk flow gra- 
dient. Note also that the z-component never crosses the 
y = line which indicates the formation of a disk-like flow 
structure with rotation in an unchanging sense. However, 
I do not expect the specific angular momentum to reach 
a quasi-steady state because of the three-dimensional na- 
ture of my simulations: angular momentum can continue 
to be accreted from directions outside the plane of the disk 
and thus probably disturbing the disk, too. 



4. Comparison of accretion rates 

4-1- Mass accretion rates 

I collected all the means of the mass accretion rates and 
their standard deviation into Fig. [l^ together with the ap- 
propriate values for models without gradients. The stan- 
dard deviation is shown as "error bar" in order to give an 
indication on how precisely one should interpret the mean 
values. 

(a) The mass accretion rates are fairly independent of 
the density gradient strength for small 7, but a decrease of 
the rates might be present when increasing the gradients: 
a slight trend toward smaller rates for larger gradients 
seems possible, (b) A clear increase of the rates is visible 
when decreasing the index 7. When going from 5/3 to 
1.01 the accretion rate increases by over a factor of two. 
(c) Models with smaller Mach numbers have larger rates. 
This trend also applies to models with a velocity gradient 
as can be seen in Fig. 12 of Rl. 

Because the fluctuations of the mass accretion rate 
slightly decrease with decreasing 7, while at the same time 
the rate itself increases, it follows that the relative fluctu- 
ation decreases strongly with decreasing 7 as can be seen 
in Fig. 14l No obvious trend seems visible on how the rel- 



increase might be present when comparing the larger gra- 
dient fluctuation with the smaller gradient ones. 

The dependence of the relative fluctuations on velocity 
seems to invert: while they are clearly larger for smaller 
Mach numbers for the 7 = 1.01 models, they are smaller 
for smaller Mach numbers in the 7 = 5/3 models. And 
last, the simulations with 7 = 4/3 produce relative fluc- 
tuations that vary more strongly with gradient strength 
than with Mach number. 

4-2. Specific angular momentum 

As has already been described in Sect. ||the specific angu- 
lar momentum reaches the analytic values given by Eq. (0) 
to within 10%, but only as long as the accretion flow is 
roughly stable and only for the models with small density 
gradients, e p = 0.03. In Fig. [l5] I show the specific angu- 
lar momentum averaged over the whole time the models 
were evolved, excluding the initial transients. This gives 
a better picture of how the very unstable flow tends to 
decrease the average and produce a large fluctuation of 
the specific angular momentum. So I included as "error 
bars" the standard deviation of the fluctuations around 
the mean. 

The unit I chose in Fig. [l^ is the specific angular mo- 
mentum that a vortex just at the surface of the accretor 
would have if it spun with the local Kepler velocity. Thus 
unity in these units indicates a Kepler orbit and mat- 
ter that has more specific angular momentum than this 
would be flung off. One can see in Fig. |l5| that the matter 
accreted is well below this value. For comparison, I also 
plot, using plus signs (+), the specific angular momen- 
tum as given by Eq. (0). These values tend to lie above 
the mean, especially for the models with large gradients 
(N,Q,U), and are well within one standard deviation from 
the mean for the small gradient models (M,P,T). For the 
latter models the fluctuations are so large that occasion- 
ally the sign of the accreted specific angular momentum 
reverses; this is indicated by the "error bars" extending 
completely down to the x-axis. 

If one is only interested in the average specific angu- 
lar momentum that is accreted, an inspection of Table [l] 
yields that the whole range of values between zero and 
about 70% is attained depending on the model parame- 
ters. Livio et al. (1986) reported values between 10% and 
20% for their parameters. 

If one reduces the size of the accretor, some point will 
be reached when the maximum specific angular momen- 
tum that can be accreted will become smaller than the 
amount present in the accretion cylinder. Setting Eq. (0) 
equal to Eq.(10) in Rl a relation is obtained between the 
gradients in the flow (e) and the radius of the accretor R±: 



ative fluctuations change with gradient strength; a slight ~ ^ ^ £v £p ^ ^ A 
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So for accretors smaller than this radius the angular mo- 
mentum accretion should no longer be dominated by what 
is given in the accretion cylinder. For e v =0 and e p =0.03 
and 0.2 we obtain R* « 3 ■ KT 5 i? A and « 10" 3 i? A . 
Both these values are below what is currently possible 
to simulate numerically, but could be important in astro- 
physical objects. On the other hand for model VS, which 
has e v =0 and e„=1.0, an accretor radius of 0.03i?A re- 
sults, which is larger than the numerically used radius 
of 0.02i?A- Thus in model VS not all the specific angular 
momentum offered in the incoming bulk flow can be ac- 
creted. Vice-versa, the maximum gradient e that can be 
accommodated by an accretor of given radius R± is 

_max / <3o .max _ /16 R+ . , 

The right panel of Fig. ^ compares the specific angular 
momentum accreted between the density gradient models 
presented here and the velocity gradient models shown in 
Rl. The ratio (cf. legend of the x-axis of Fig. ||) is plot- 
ted of the specific angular momentum (values given in 
Table |l|) for these two sets of models. If a pair of mod- 
els experiences the same reduction in accretion of specific 
angular momentum, the ratio plotted would be one. If 
the density-gradient model of the pair suffers a greater 
decrease of accretion (due to e.g. stronger relative fluctu- 
ations of the unstable flow) this will be reflected by a ratio 
that is smaller than unity. 

The pair MF/JS is less than zero, because the aver- 
age angular momentum accreted by model MF is actually 
retrograde to the bulk flow momentum. For this pair, as 
well as for MS/IS, the very large 'error' bar indicates that 
the fluctuations due to the unstable flow are very large 
compared to the average and so the latter value does not 
permit a strong statement. So although for four of the five 
models the ratios are above unity, which would indicate 
that the unstable flow actually increases the angular mo- 
mentum accretion for the density gradient models as com- 
pared to the velocity gradient models, this reasoning is not 
a credible one. Additionally one has to keep in mind that 
the specific angular momentum of the incoming bulk flow 
of models KS and LS is actually larger than the maximum 
permitted by the Eq. [ll]. Thus the accreted momentum 
will be decreased due to this effect, too (angular momen- 
tum barrier), which explains why these model-pairs have 
the largest ratios. 

4-3. Correlations 

A correlation was found in Rl (Fig. 17) between the mass 
accretion rate and the specific angular momentum: the 
rate decreases when the magnitude of the momentum is 
largest. In Fig. |l6| I draw a similar plot as in Rl, each dot 
connecting the two quantities for every second timestep of 
the numerical simulation. A similar trend can be seen for 
model QS, which confirms that the dynamics is similar: 



when the flow does not rotate (low specific angular mo- 
mentum) if falls down the potential to the surface of the 
accretor and thus produces a higher mass accretion rate. 
For all other models a correlation is not obvious. 

5. Summary 

For the first time a comprehensive numerical three- 
dimensional study is presented of wind- accretion with a 
density gradient using a high resolution hydrodynamic 
code. I vary the following parameters: Mach number of 
the relative flow (Mach AA=3 and 10), strength of the 
density gradient perpendicular to this flow (e p =3%, 20% 
and 100% over one accretion radius), and adiabatic in- 
dex (7=5/3, 4/3, and 1.01). The results are compared (a) 
among the models with different parameters, (b) to some 
previously published simulations of models with density 
and velocity gradients, and (c) also to the analytic esti- 
mates of the specific angular momentum. 

1. All models exhibit active unstable phases, which are 
very similar to the models without gradients. Only the 
mildly supersonic case (M.=1A) displayed a relatively 
steady flow (as compared to the faster flow models). 
The accretion rates of mass, linear and angular mo- 
mentum fluctuate with time, although not as strongly 
as published previously for 2D models. 

2. Depending on the model parameters, the average spe- 
cific angular momentum accreted is roughly between 
zero and 70% of the analytical estimate, which assumes 
that all angular momentum within the accretion cylin- 
der is actually accreted. 

3. The mass accretion rates of all models with density 
gradients are equal, to within the fluctuation ampli- 
tudes, to the rates of the models without gradients 
(published previously), although the accretion rates 
might seem to decrease slightly when increasing the 
density gradient. The fluctuations of the mass ac- 
cretion rate in all models hardly vary with gradient 
strength. 

4. The overall qualitative flow dynamics as well as the 
mass accretion rates are very similar to what has been 
published on models with velocity gradients. Of course 
the accretion rate of angular momentum and its spe- 
cific values (i.e per mass unit) are reduced if one com- 
pares equal gradients for both velocity and density, 
well in accordance with the analytic estimates. This re- 
duction means that in the density gradient models the 
fluctuations due to the unstable accretion flow have a 
greater influence on the angular momentum accretion 
than for the velocity gradient models. 

5. The models with small gradients (e p =0.03) display an 
initially quiet stable phase, in which the specific an- 
gular momentum of the matter accreted is within 10% 
of the analytic estimate. Thus for the quiescent phases 
the analytic values are appropriate. The average drops 
when the flow becomes unstable. 



18 



M. Ruffert: Non-axisymmetric wind-accretion simulations 



6. The model with very large density gradient (e p = 100% 
over one accretion radius) was the only one for which 
the accreted angular momentum was always prograde 
with respect to the angular momentum available in the 
incoming flow. Here the amplitude of the perturbation 
due to the unstable flow is much smaller than the aver- 
age angular momentum accreted. However the specific 
angular momentum of the incoming flow in this case is 
larger than the maximum given by the Kepler velocity 
times the radius of the accretor surface. 
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